Integrating atomistic molecular dynamics simulations, experiments, and network analysis to study protein dynamics:strength in unity

In the last years, we have been observing remarkable improvements in the field of protein dynamics. Indeed, we can now study protein dynamics in atomistic details over several timescales with a rich portfolio of experimental and computational techniques. On one side, this provides us with the possibility to validate simulation methods and physical models against a broad range of experimental observables. On the other side, it also allows a complementary and comprehensive view on protein structure and dynamics. What is needed now is a better understanding of the link between the dynamic properties that we observe and the functional properties of these important cellular machines. To make progresses in this direction, we need to improve the physical models used to describe proteins and solvent in molecular dynamics, as well as to strengthen the integration of experiments and simulations to overcome their own limitations. Moreover, now that we have the means to study protein dynamics in great details, we need new tools to understand the information embedded in the protein ensembles and in their dynamic signature. With this aim in mind, we should enrich the current tools for analysis of biomolecular simulations with attention to the effects that can be propagated over long distances and are often associated to important biological functions. In this context, approaches inspired by network analysis can make an important contribution to the analysis of molecular dynamics simulations

CAncer bioMarker Prediction Pipeline (CAMPP) - A standardized framework for the analysis of quantitative biological data

With the improvement of -omics and next-generation sequencing (NGS) methodologies, along with the lowered cost of generating these types of data, the analysis of high-throughput biological data has become standard both for forming and testing biomedical hypotheses. Our knowledge of how to normalize datasets to remove latent undesirable variances has grown extensively, making for standardized data that are easily compared between studies. Here we present the CAncer bioMarker Prediction Pipeline (CAMPP), an open-source R-based wrapper (https://github.com/ELELAB/CAncer-bioMarker-Prediction-Pipeline -CAMPP) intended to aid bioinformatic software-users with data analyses. CAMPP is called from a terminal command line and is supported by a user-friendly manual. The pipeline may be run on a local computer and requires little or no knowledge of programming. To avoid issues relating to R-package updates, a renv .lock file is provided to ensure R-package stability. Data-management includes missing value imputation, data normalization, and distributional checks. CAMPP performs (I) k-means clustering, (II) differential expression/abundance analysis, (III) elastic-net regression, (IV) correlation and co-expression network analyses, (V) survival analysis, and (VI) protein-protein/miRNA-gene interaction networks. The pipeline returns tabular files and graphical representations of the results. We hope that CAMPP will assist in streamlining bioinformatic analysis of quantitative biological data, whilst ensuring an appropriate bio-statistical framework

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally. DOI: http://dx.doi.org/10.7554/eLife.17505.00

Inhibitors of the Cdc34 acidic loop:a computational investigation integrating molecular dynamics, virtual screening and docking approaches

AbstractAmong the different classes of enzymes involved in the ubiquitin pathway, E2 ubiquitin-conjugating enzymes occupy a central role in the ubiquitination cascade. Cdc34-like E2 enzymes are characterized by a 12–14 residue insertion in the proximity of the catalytic site, known as the acidic loop. Cdc34 ubiquitin-charging activity is regulated by CK2-dependent phosphorylation and the regulatory mechanism involves the acidic loop. Indeed, the phosphorylation stabilizes the loop in an open conformation that is competent for ubiquitin charging.Cdc34 is associated with a variety of diseases, such as hepatocellular carcinomas and prostatic adenocarcinomas. In light of its role, the discovery of potential inhibitory compounds would provide the mean to effectively modulate its activity.Here, we carried out a computational study based on molecular dynamics, virtual screening and docking to identify potential inhibitory compounds of Cdc34, modulating the acidic loop conformation. The molecules identified in this study have been designed to act as molecular hinges that can bind the acidic loop in its closed conformation, thus inhibiting the Cdc34-mediated ubiquitination cascade at the ubiquitin-charging step. In particular, we proposed a pharmacophore model featuring two amino groups in the central part of the model and two lateral aromatic chains, which respectively establish electrostatic interactions with the acidic loop (Asp 108 and Glu 109) and a hydrogen bond with Ser 139, which is one of the key residues for Cdc34 activity

Computational Structural Biology of S-nitrosylation of Cancer Targets

Nitric oxide (NO) plays an essential role in redox signaling in normal and pathological cellular conditions. In particular, it is well known to react in vivo with cysteines by the so-called S-nitrosylation reaction. S-nitrosylation is a selective and reversible post-translational modification that exerts a myriad of different effects, such as the modulation of protein conformation, activity, stability, and biological interaction networks. We have appreciated, over the last years, the role of S-nitrosylation in normal and disease conditions. In this context, structural and computational studies can help to dissect the complex and multifaceted role of this redox post-translational modification. In this review article, we summarized the current state-of-the-art on the mechanism of S-nitrosylation, along with the structural and computational studies that have helped to unveil its effects and biological roles. We also discussed the need to move new steps forward especially in the direction of employing computational structural biology to address the molecular and atomistic details of S-nitrosylation. Indeed, this redox modification has been so far an underappreciated redox post-translational modification by the computational biochemistry community. In our review, we primarily focus on S-nitrosylated proteins that are attractive cancer targets due to the emerging relevance of this redox modification in a cancer setting

ENCORE:Software for Quantitative Ensemble Comparison

There is increasing evidence that protein dynamics and conformational changes can play an important role in modulating biological function. As a result, experimental and computational methods are being developed, often synergistically, to study the dynamical heterogeneity of a protein or other macromolecules in solution. Thus, methods such as molecular dynamics simulations or ensemble refinement approaches have provided conformational ensembles that can be used to understand protein function and biophysics. These developments have in turn created a need for algorithms and software that can be used to compare structural ensembles in the same way as the root-mean-square-deviation is often used to compare static structures. Although a few such approaches have been proposed, these can be difficult to implement efficiently, hindering a broader applications and further developments. Here, we present an easily accessible software toolkit, called ENCORE, which can be used to compare conformational ensembles generated either from simulations alone or synergistically with experiments. ENCORE implements three previously described methods for ensemble comparison, that each can be used to quantify the similarity between conformational ensembles by estimating the overlap between the probability distributions that underlie them. We demonstrate the kinds of insights that can be obtained by providing examples of three typical use-cases: comparing ensembles generated with different molecular force fields, assessing convergence in molecular simulations, and calculating differences and similarities in structural ensembles refined with various sources of experimental data. We also demonstrate efficient computational scaling for typical analyses, and robustness against both the size and sampling of the ensembles. ENCORE is freely available and extendable, integrates with the established MDAnalysis software package, reads ensemble data in many common formats, and can work with large trajectory files